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Abstract: We study a parton shower description, based on a dipole picture, of the fi- 
nal state in electron-positron annihilation. In such a shower, the distribution function 
describing the inclusive probability to find a quark with a given energy depends on the 
shower evolution time. Starting from the exclusive evolution equation for the shower, we 
derive an equation for the evolution of the inclusive quark energy distribution in the limit 
of strong ordering in shower evolution time of the successive parton splittings. We find 
that, as expected, this is the DGLAP equation. This paper is a response to a recent paper 
of Dokshitzer and Marchesini that raised troubling issues about whether a dipole based 
shower could give the DGLAP equation for the quark energy distribution. 
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1. Introduction 

Parton shower algorithms, when coupled to hadronization models, provide a way of gener- 
ating cross sections for exclusive final states according to approximations based on QCD. 
One can use these algorithms also to produce predictions for inclusive observables. Of 
special interest are predictions that, in a perturbative expansion, involve large logarithms 
of ratios of different momentum scales in the physical problem. In many cases, there are 
also predictions based on the full field theory. In these cases, it is of interest to know if a 
particular parton shower algorithm gives results consistent with the full field theory. 

In this paper, we investigate one such case. We consider the inclusive parton energy 
distribution in electron-positron annihilation: the probability that a parton carries a frac- 
tion x of the maximum possible energy of a parton. (We define the energy distribution 
precisely in Sec. 2.) We ask how the parton energy distribution changes as the resolution 
with which we look at the shower changes. We use the shower evolution equation of Ref. [1] , 
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which may be said to describe a dipole shower in that the emission of a gluon from a parton 
I can interfere with the emission of a gluon from another parton k. The shower evolution 
equation describes the evolution of exclusive parton states. We manipulate the equation 
to produce an evolution equation for the inclusive parton energy distribution. This energy 
distribution of the partons should obey, approximately, the lowest order DGLAP evolution 
equation [2] (like parton distribution functions). 1 

This investigation is prompted by a recent paper of Dokshitzer and Marchesini [4], 
which raises some troubling issues with respect to the flow of momentum in the parton 
splittings. Ref. [4] suggests that perhaps a dipole shower cannot get the evolution of the 
energy distribution right. This is a significant issue since parton shower Monte Carlo event 
generators are of great importance in the design and interpretation of experiments and 
because some of these are based on a dipole picture. For instance, the program Ariadne 
[5] is based on dipoles and, in fact, introduced the dipole idea as an organizing principle for 
parton showers. 2 In particular, the recent program of Dinsdale, Ternick and Weinzierl [6] 
and the program of Schumann and Krauss [7] are both of the dipole type. These programs 
are based on the Catani-Seymour dipole splitting functions and momentum mapping [8] 
that are often used in next-to-leading order perturbative calculations and can be applied to 
showers as described in Ref. [9]. Furthermore, the final state shower in version 8 of Pythia 
[10, 11] is effectively a Catani-Seymour dipole shower. Thus if dipole showers cannot work 
for as significant a quantity as the evolution of the parton energy distribution, we ought to 
know about it. 

For our investigation, we consider the process e + + e~ — > partons in a shower picture. 
Our principle example is the shower evolution equation of Ref. [1] . The organizing principle 
of the shower is that one starts at the hard interaction 7 — ► q + q and moves to softer 
interactions, always factoring the softer interaction from previous harder interactions. We 
find in Sec. 3 and Appendix A that with the definitions of Ref. [1], we do obtain the 
expected DGLAP equations. We will also show in Sec. 5 that, at least for the case of 
q — > q + g splittings, the DGLAP equation results if one uses the splitting functions and 
momentum mapping of the Catani-Seymour dipole scheme. For either choice of splitting 
functions and momentum mapping, we take virtuality as the measure of hardness. One 
could also use the transverse momentum of the splittings to order successive splittings as 
in Refs. [6, 7, 10, 11]. As we outline in Sees. 3.6 and 5, transverse momentum ordering 
also produces the DGLAP equation. We devote Sec. 4 to the role of soft gluons in the 
evolution equation, since this was the point of interest in Ref. [4] . 

The main point of this paper is to show that one can obtain the lowest order DGLAP 
equation from the shower evolution equation of Ref. [1] in a simple fashion, using the essen- 
tial approximation that the shower splittings, which are by definition ordered in virtuality, 
are strongly ordered. In fact, successive splittings are strongly ordered if the splitting 
probability for each step dq 2 in virtuality is small, which follows if a s is small. 

Before beginning the main analysis, we provide an informal and qualitative description 
of what we think the issue raised in Ref. [4] is. Suppose that in the first step of a shower 

lr rhe application to parton decay functions can be found in Ref. [3]. 

2 However, we do not investigate whether the DGLAP equation results in Ariadne. 
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Figure 1: Emission of a soft gluon from a nearly collinear quark-gluon jet. The gluon is at a wider 
angle to the jet than the jet opening angle, but is very soft. The gluon does not resolve the jet, 
so the emission is as if the soft gluon were emitted from an on-shell massless parton moving in the 
direction of the jet and carrying the net color of the jet, as indicated on the right hand side of the 
approximation. 

description for e + + e~ —> partons, the initial quark emits a gluon with an energy u)\ equal 
to 1/3 of the initial quark energy, E q . Then the daughter quark has energy 2/3 x E q . 
Suppose that the emission angle 9\ is very small, say 9\ = 0.01. Now suppose that a soft 
gluon with energy loi is emitted from the daughter quark or the daughter gluon at an angle 
02 that is large compared to 9\, say 62 = 0.1. This is the second emission in a virtuality 
ordered shower, so the virtuality of the emission, q\ ~ E q u)2Q\ must be smaller than the 
previous virtuality, q\ = E q iO\Q\. Thus we need ix>2 < w\Q\jQ\ ~ 10~ 2 wi. 




Figure 2: Emission of a soft gluon from a nearly collinear quark-gluon jet when the gluon in the 
jet is considered as 3 ® 3 state. Emission from the narrow dipole cancels. 

Now we can ask where the energy of the soft gluon comes from. We note that the 
soft gluon cannot resolve the narrow jet. Thus, as indicated in Fig. 1, the matrix element 
for emitting the soft gluon is approximately the same as the matrix element for emitting 
the soft gluon from an on-shell quark moving in the direction of the mother quark. In an 
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angular ordered shower, the emission is, in fact, modeled in this way [12]. This suggests 
that the energy of the soft gluon should come from the mother quark, so that 2/3 of it is 
taken from the daughter quark and 1/3 from the daughter gluon. However, this suggestion 
is (at least within the intuitive picture used here) misleading. If we make a simple leading 
color approximation, then the daughter gluon carries color 3 <8> 3 instead of color 8. The 
soft gluon could be emitted from the 3 color in the daughter gluon, but the soft gluon 
cannot resolve this color from the 3 color of the daughter quark. Thus this part of the 
gluon emission approximately cancels, as depicted in Fig. 2. Not all of the soft gluon 
emission amplitude cancels. The amplitude for emitting the soft gluon from the 3 part of 
the daughter gluon, as depicted in Fig. 3, remains. This amplitude is approximately the 
same as that for emitting the soft gluon from the mother quark, as in the right hand side 
of Fig. 1. However, the soft gluon emission comes after the mother quark has split into a 
daughter quark and a daughter gluon. We thus see that it is the daughter gluon that gives 
up part of its energy to supply the energy of the soft gluon. 

Does this issue matter? Perhaps, but it does not matter to the evolution of the energy 
distribution of the quark. The reason is that the energy taken by the soft, wide angle gluon 
is too tiny to make a difference in the limit in which we can regard the virtualities in the 
shower as strongly ordered. We will see this point appear in the actual shower evolution 
equations in the following sections. 




Figure 3: Emission of a soft gluon from a nearly collincar quark-gluon jet when the gluon in the 
get is considered as 3 ® 3 state. Emission from the part of the gluon that makes part of a wide 
angle dipolc remains. 

We point out that the virtuality ordering matters to this argument. Consider the case 
that u>2 < wi but L02 > LO\0\jQ\. Then the soft gluon emission is harder than the collinear 
q — > q + g splitting, so that the soft gluon emission belongs first in the shower, and part of 
its energy is to be deducted from what the daughter quark would otherwise get. Suppose 
that one were to use an energy ordered shower, despite the fact that the approximations 
used in factoring the softest splitting from harder splittings require that softer splittings 
are later in the shower. Then in a dipole style shower this soft gluon emission would come 
second and the non-cancelling contribution would come from Fig. 3. The wrong parton 
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would then be giving up its energy to supply the soft gluon energy. One would then not be 
surprised to not get the DGLAP evolution equation from such a shower evolution if one did 
not neglect u>2 compared to E q . We suspect that something like this may have happened 
in Ref. [4]. In this paper, we use a virtuality ordered shower, as in Ref. [1]. 

2. The parton energy distribution 

We consider the process e + + e~ — > partons in a shower picture, using the notation of 
Ref. [1]. The shower evolution equation of Ref. [1] includes quantum interference among 
colors and spins. We note that this shower evolution equation is not directly practical for 
generating events. In order to implement the shower evolution equation, one would need 
an approximation scheme such that the approximated shower evolution can be practically 
implemented and such that the approximations can be systematically improved at a cost 
of computer time. In the lowest order approximation, one averages over spins and takes 
the leading color approximation, which yields an evolution equation that can be written 
as a Markov process [13]. We will discover in this paper that we get the conventional 
DGLAP evolution equation for the energy distribution of partons as a function of resolution 
either with the full treatment of color and spin or in the leading color, spin averaged 
approximation. 

In the notation of Ref. [1] , states in the sense of statistical mechanics are represented 
by ket vectors |p), while possible measurements are represented by bra vectors (i* 1 . Thus 
(-F|/?) is the probability that one obtains a particular result F from a measurement on an en- 
semble of systems represented by |p) . We use basis states labelled by lists {p, f, s', c', s, c} m 
of parton quantum numbers. Here {pi, . . . ,p m } = {p} m are the momenta of the partons 
and {f} m are their flavors. The total momentum is Q = YlPi- The spins are specified by 
{s 1 , s} m , with two spin labels for each parton because we use the quantum density matrix 
in spin space in order to represent possible interference among spin states. The colors 
are similarly specified by {c',c} m . Quantum spin and color states are represented by bra 
and ket vectors with angle brackets, as in the quantum spin inner product ({s'} m |{s} m ). 
A much more detailed specification of our notation is provided in Ref. [1]. In Ref. [1], 
the quarks can have non-zero masses, but in this paper we take all of the partons to be 
massless. 

Let \p(t)) = U{t, 0) p(0)) be the statistical state at shower time t. Here t = at the 
beginning of the shower, which starts with the hard process, and t — > oo gives the statistical 
state after hadronization. The total cross section is 

( 1 \p( t )) = J2^ / [d{ P ,f,s',c',s,c} m ] 

x (1|{P, /, s',c',s,c} m )({p,f,s',c',s,c} m \p(t)) 

1 f ' ' 

= E rn T / HP'/' 5 '' '' 5 ' }™] ({ s '}m\{s}m) ({c'} m \{c} m ) 

X {{P,f,s',C,S,c} m \p(t)) . 
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This is independent of t because (l \U(t, 0) = (l|. We now introduce a measurement 
function (x,/ a | that measures the probability to find in the final state a parton (or a 
hadron in the case t — > oo) with flavor / a having a fraction x of the maximum possible 
energy for a single particle, \J Q 2 /2. The definition is 

{x,U\p(t))= j- / [d{p,/y,</,s,c} m ] 

X (s, /a|{P, /, s', c', S, C} m ) ({p, /, s', c', S, c} m |p) 

= f / [ d {P 5 /> s ' 1 c/ > s ) c }m]<{s / }m|{s}m} <{c'} m |{c} m ) ( 2-2 ) 

m ' ^ 

X ^2 S U,fi S(x- 2P Q2^ {{p,f,s',C,S,c} m \p(t)) . 

The distribution measures all partons, so there is a sum over parton labels i. The energy 
fraction x for parton i is defined to be 2pi • Q/Q 2 . Note that if we integrate over x with a 
factor x and sum over flavors, we get twice the total cross section: 

^2 I dxx (x,f a \p(t)) = / KP./.S'.C'-S.C}™] ({s'}m\{s}m) ({c'} m \{c} m ) 



X ^2^ ^ ({P)/) s/ > c/ ) s > c }m|p) 

= 2(l|p(t)) . 

(2.3) 

When we apply this at t = oo, we have a physical observable giving, say, the energy 
distribution of hadrons of flavor /h in an electron-positron collision, (x, /h|/o(oo)). This 
observable can be written in terms of the cross section to produce a parton of flavor / a 
with a momentum fraction x at a resolution scale t convoluted with the parton decay 
function that gives the probability for that parton to decay to a hadron of the desired 
flavor carrying a fraction x/x of the parton's momentum: 

(x,/ h | P (oo)) « d i D kiu{^ t ) MM) ■ ( 2 - 4 ) 

/a 

This equation is approximate because factoring out the parton decay function involves 
neglecting power suppressed contributions. The parton decay functions depend on the 
resolution scale t that we used to factor the cross section. This dependence is given by 
the DGLAP equation [3], Since the observable (x, /h p(oo)) does not depend on t, the 
partonic function (x, / a |/o(t)) must also depend on t and this dependence must be given 
approximately by the DGLAP equation. We will investigate this evolution scale dependence 
within our dipole-like equations for shower evolution. 

3. Shower evolution for the energy distribution 

We seek to discover how (x,/ a |/?(t)) evolves with shower resolution parameter t. We will 
consider the case in which / a is a quark flavor, / a G Q = {it, d,s, . . .}. This was the case 
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examined in Ref. [4]. The case of an antiquark is essentially identical. The case of a gluon 
is analyzed in Appendix A. 

From the evolution equation [1] for p(t)), we have 



j t (x,U\p(t)) = (x,U\Hi(t)-V(t)\p(t)) . 



(3.1) 



Here 7~Ci(t) is the parton splitting operator and V(t) is the virtual "no-splitting" operator 
that insures that the probability for a splitting plus the probability for no splitting equals 
unity. We want to investigate [x, f a \7ii(t) — V(t) in the limit of strongly ordered 

virtualities in the shower. That is, the virtuality in the shower splitting described by 
TCl(t) — V(t) should be small compared to the virtualities in the state \p(t)) before the 
splitting. We will make approximations that reflect this condition. 

When we insert a complete sum over statistical states next to \p(t)) in Eq. (3.1), we 

get 

^(£>/a|p(*)) = f [ d {PJi s '> c 'i S ' C }rn] 

m 

x (x, f a \Hi(t) - V(t)\{p, /, s', c', 8, c} m ) 
x {{P,f,s',c',s,c} m \p(t)) . 

In the following subsections, we will evaluate the matrix element of Tii(t) — V(t) in the 
approximation of strongly ordered virtualities. In the last subsection, we will insert the 
limiting result for this matrix element back into Eq. (3.2). This will yield the DGLAP 
equation. 

In the analysis of this section, we will freely make approximations that follow from 
having strongly ordered virtualities. Some of these approximations, related to soft, wide 
angle gluon emissions, were the main focus of Ref. [4]. In Sec. 4, we turn to a detailed 
examination of these soft emission approximations. 

3.1 The splitting matrix element 

To see what Eq. (3.2) means, we examine the matrix element of 7i\{t) — V(t) using the 
results in Sec. 12 of Ref. [1], 

(x, f a \Hl(t) - V(t)\{p,f,s',C,S,c} m ) = 

dC P 0(( p e m P } m , Cf)) Sft - 7}({p, f} m+l ) 



i he® i(fi) 

■ m+l 



X 



2pj • Q 

Q 2 



2 Pi -Q 

Q 2 



X lo(fm+l T R ({c'} m \{c} m ) ({s'} m \{s} m ) Wa({f,p} m+1 ) 

+ 6(f m +l = g) E ({ c '}™\tk(fk ~+ fk + /m+l) 4(/z Si + L+l)\{ch 
k+l 



(3.3) 



x ({s'} m \{s} r 



2Ai k ({p} m+1 )wi k ({f,p} m +l) - W U ({f,p}m+l) 
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This formula is a little complicated because it contains several ingredients. However, its 
important features are easily understood. 

In the second line, there is a sum over the index I of the parton that splits. The 
daughter partons are labelled I and m + 1. Then there is a sum over the flavors of the 
daughter partons into which parton I splits. These flavors are denoted by Cf = (fi, fm+l)- If 
fi is a quark or antiquark flavor, it must split into the same flavor quark or antiquark plus a 
gluon. If fi is a gluon, it can split into two gluons or a quark-antiquark pair. Then there is 
an integration over splitting variables that determine the momenta of the daughters. These 
can be a virtuality y, a momentum fraction z, and an azimuthal angle <j). This integration 
is indicated as d£ p . The momenta pi of the partons after the splitting are given in terms 
of the momenta pi before the splitting together with the splitting variables £ p . Finally in 
the second line there is a delta function that defines the shower time t, for which we will 
use the virtuality of the splitting, 

Ti({pJ} m +i)= log (— ) . (3.4) 

\2pi -p m+ ij 

In the third line there is a difference of two terms. The first, with /j and pi, comes 
from TLiit). The second, with /j and pi, comes from V(t). In this virtual contribution, 
the energy fraction x is set to the energy fraction of parton i before the splitting and we 
integrate inclusively over the splitting variables. 

The functions that follow are defined in Ref. [1]. They are derived from the spin 
dependent splitting functions for the splitting of a parton with label I to produce two 
daughter partons with, in our labeling scheme, labels I and m + 1. 

The line starting 0(/ m +i / g) is for a g — > q + q splitting. It has a color factor 
^R({ c '}m|{ c }m) with Tr = 1/2. The function wu({f,p} m+ i) is the spin-averaged split- 
ting function for g — ► q + q, as given in Ref. [1]. This term is depicted in Fig. 4. The 
dependence on the mother parton spin is trivial, ({s'} m |{s} m ), because we sum over the 
daughter parton spins and integrate over the azimuthal angle of the splitting, as indicated 
in Eq. (12.18) of Ref. [1]. 

The lines starting 6{f m +i = g) are for a q — > q + g splitting. It has a nontrivial color 
structure ({c'} m |ifc t||{c} m ). This corresponds to an interference graph and is the color 
factor for emitting a gluon from parton I in the quantum amplitude and absorbing it on line 
k in the complex conjugate amplitude. 3 This color factor multiplies a sum of two terms. 
The first, containing wik, corresponds to the interference graphs, depicted in Fig. 5. The 
second, containing wu, corresponds to the direct graphs, depicted in Fig. 6. The functions 
wik and wu are given in Ref. [1], while function A^, or rather possible choices for it, are 
given in Ref. [14]. 

3.2 Strongly ordered virtualities 

We can immediately simplify Eq. (3.2) in the strong ordering approximation. There are 
sums over partons i that participate in the measurement represented by (x,/ a |. Thus x 

3 This factor is written as the matrix element of ■ Tj in the work of Catani and Seymour [8]. 
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Figure 4: A graph contributing to wn for a g — > g + g splitting. 




Figure 5: An interference graph contributing to wik for aij^g + g splitting. 



is set equal to 2pi ■ Q/Q 2 in the term from Hi(t) and 2pi ■ Q/Q 2 in the term from V(t). 
Consider the contribution if i is the same in both terms and is neither of the indices of the 
daughter partons, iorm+1. For a splitting that is nearly collinear or is very soft, or both, 
we have pi ~ pi for such a spectator parton. In fact, we would have pi = pi except for the 
need to take some momentum from the spectator partons in order to balance momentum in 
the splitting while keeping both the daughter partons and the mother parton on shell. The 
amount of momentum needed tends to zero as the virtuality in the splitting approaches 
zero. For this reason, pi — pi —* in the limit of small virtuality in the splitting. Thus in 
the limit of strongly ordered virtualities in the shower, we drop the terms for i not equal 
to I or m + 1. This becomes yet simpler because we consider the case that f a is a quark 
flavor. With our labeling scheme, parton m + 1 is always either a gluon or (for g — ► q + q) 
an antiquark. Thus, in fact, the only term that remains in the sum over i is that for i = I. 
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This gives 



[x, U\Hi{t) - V(t)\{p, f, s', c', s, c} m ) « 

Yl j d ^ *«P e r l(iP}m, Cf)) S(t - r,({p, /}m+l)) ({4m|{4m) 



i Cf6*l(/i)' 



5; f (5| I 

JSLiJl 



2pi-Q 
Q 2 



St r 5\ x 



2pi-Q 

Q 2 



X { 0(fm+l + g) T R ({c'}m\{c} m ) Wu({f,p}m+l) 



(3.5) 



+ 0(fm+l = g) Y ({c}m\tkUk -> fk + fm+l) t\(fl -> fl + fm+l)\{c} m ) 
k+l 

2A lk ({p} m+1 )wik({f,p} m +i) ~ Wll({f,P}m+l) 



The bookkeeping for the flavor structure can now be made more explicit. Consider the 
two terms on the third line of Eq. (3.5), representing a real splitting and a virtual splitting. 
If the mother parton I is a gluon, the virtual splitting contribution vanishes because f a ^ fi- 
In the real splitting contribution, the Kroneker delta restricts the sum over flavors Q in 
the splitting to just one term, (//, / m +i) = (/ a , —fa)- If parton I is not a gluon, it must in 
fact be a quark with flavor fi = f a . Then again, there is only one term in the sum over 
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flavors Cf in the splitting, namely (fi, f m +i) = (/a,g)- Thus 



{xJ^Hitt) - V(t)\{p,f,s',c',s,c} m ) « 

^/rfC P ^(C P erK{p} w ,Cf))^(i-7]({p,/} m+ i)) 2 <{a'} m |W, 



X i 0(fi =g,fl = /a, /m+1 = "/a 



fljrt 



x 7R<{c'} m |{c} m ) p}m+l) <sf&- 2 ^ 2 ^ 

+ = fa, fl = /a, /m+1 = g) 

X X] ({ c '}m\tk(fk -> fk + fm+l)t}(fl -+fl + /m+l)|{c} r 

2Afc({p}m+l)^fc({/,P}m+l) - %({/.PWl) 



(3.6) 



A further simplification is possible. The interference diagram of Fig. 5, represented by 
the function wik, is important in the development of the exclusive final state in a parton 
shower. However, wit is singular only in the limit of a soft gluon emission. It is not singular 
for emission with p m +i almost parallel to pi. In the limit of a soft gluon emission, the soft 
gluon does not take any momentum from parton /. Thus p\ ~ pi, so that the two terms on 
the last line of Eq. (3.6) cancel. Thus, for our inclusive observable, in the strong ordering 
limit the term in Eq. (3.6) proportional to Wik can be dropped. We recognize that this step 
touches on the main point of Ref. [4], so we will return to provide a more detailed analysis 
of it in Sec. 4. 

Once we have dropped wik, the sum over k in Eq. (3.6) is trivial. As in Ref. [8] and 
in Eq. (8.25) of Ref. [1], we use the fact that the complete quantum state is a color singlet 
to write 



^2 ({ C '}m\tk(fk ->■ fk + frn+l) t\(fl ~> fl + /m+l)|{ c }m) 

= - ({c'}m\tl(fk -> fk + frn+l)t\{fl -> fl + f m+ l)\{c} m ) 
= -C F ({c'} m \{c} m ) . 



(3.7) 



This simplification holds when we treat color exactly. It also holds when we use the leading 
color approximation, as explained in Ref. [13]. 

It will also prove useful to define x = 2pi ■ Q/Q 2 and introduce into Eq. (3.5) a factor 



/ dx s(x — ^ PI q 2 ^ 



(3.8 
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With these changes, Eq. (3.6) becomes 

(x,f a \Hl(t) -V(t)\{p,f,s',C,S,c} m ) ~ 

2 / d ^ #(C P G r,(Mm, CO) s{t - mp, /Wi)) ^ 

fit f m+1 

X ({s'}m\{s}m) ({c} m |{c} m ) ^ dx 5^X - ^2^ ) 

X =g,fl = /a, /m+1 = -/a) T K W U ({f, p} m+ l) S (x - ^~^ x 



(3.9) 



+ #(/z — /a, /i — /a, /m+1 — g 

X C F uJ/z({/,p} m+ i 



(5 ( x — — — ^ a; ) — <5(x — x) 
Pi ■ Q 



3.3 Kinematics 



At this point, we introduce the variables used in Ref. [13] to describe the splitting of parton 
I with momentum pi to produce daughter parton I with momentum p\ and daughter parton 
m + 1 with momentum p m +i- We define a virtuality variable 

2/ = ~s tt- ■ (3.10) 

We use virtuality as the evolution variable, as defined in Eq. (3.4). (In Sec. 3.6, we examine 
using transverse momentum as the evolution variable.) Accordingly, the delta function that 
sets the shower time of the splitting gives 

v = Sn e ~ l - (3 - n) 

We define a lightlike vector n\ in the pi-Q plane and use this vector to define a momentum 
fraction variable z by 

Pl • Hl (3.12) 



Pl-ni+ p m+ i ■ m 

These are the variables used in Ref. [13] except that we have interchanged z <-> (1 — z) in 
order to follow the convention for a q — > q + g splitting that z is the momentum fraction 
of the daughter quark. There is a third splitting variable, the azimuthal angle of the 
transverse part of p m +i i n a reference frame in which pi and n; define the 0-3 plane. 

We can work out from the definition in Ref. [1] that if we use these splitting variables, 

dC P 9(C p emp} m ,C{))--- =^J Q V ™dy\£dz y^-.. . (3.13) 
Here, 4 

\ = ^{l + yY-Ay/x , 

1 fl V (3-14) 



1 

x V x 



4 In this paper, we denote 2pi ■ Q/Q 2 = x. In Ref. [1] we wrote 2pi ■ Q/Q 2 = 1/a/. 
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We take note that x < 1 so that A < 1. The maximum value of y is the value that makes 
A = 0. Also 1 — A is proportional to y for y — > 0. Using the delta function that sets y in 
terms of t, we have 



j <iCp«(Cp e n({ P } m x,)) sit - T,mf} m+ i) 



y f(y < 2/max) A dz 



OK Jo J 27T 

where now ?/ is given by Eq. (3.11). When we take the limit y — > 0, we have 0(y < y m ax) - > 
1. We also have A — > 1. 

The delta function 6(x — (pi • Q/pi • Q) x) relates x, = 2pi • Q/Q 2 to x = 2pi • Q/Q 2 
and the splitting variables y and z according to the momentum mapping in Ref. [1]. The 
relation is determined from 

frQ . + l v -(,-l)ll-X). (3.16) 



PlQ 2 

We will use the full relation in Sec. 4. In this section, we immediately take the limit of 
strongly ordered virtualities, y — ► 0, 



Pi • Q 
Pi ■ Q 

Thus 

{x, h\Hi(t) - V(t)\{p,f,s',c,s,c} rn ) a 
sr^Pi-Q f 1 , fd(f> sr 

SlJm + l 



z . (3.17) 



2pi-Q 



x (W}m\{s}m) ({c'} m \{c}m) J dx 5 (^X 

X yifl =g,fl = fa, fm+1 = - fa) T R W U ({f,p} m+ i) - Si X 



Q 2 J (3.18) 



+ Q(fl — /a, fl — /a, fm+1 — g) 

x C F wu({f,p} m+ x) 



- S[ x ) — S(x — x) 

z \ z 



3.4 Splitting functions 



We can take the functions wu from Ref. [13]. They have the form of a simple prefactor 
times a function of the splitting variables y and z and the mother parton energy fraction 

x, 

_ , r i \ 4vra s (F q/ (y,z,x), (f u f m+1 ) £ {(u, u), (d, d), . . . } 

^(l f 5 Prm+i) = ^ — x< * * • (3.19) 

Pl-Qy \F q/q (y,z,x), (/,,/ m+1 )e{(u,g),(d,g),...} 
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The functions F are, from Eqs. (A.l) and (2.23) of Ref. [13], 
F q / g {y,z,x) = z 2 + (1 - z) 2 , 



F q/q (y,z,x) 



(1 + X + yl 
4 A 



25 



2y 



1-5 x(l - 5) 2 (1 + y) 2 



(3.20) 



+ -(l-z)[l + y + A] , 



where 5 



z(y,z,x) 



\z + ±(l + y-\) 
1 + 1/ 



(3.21] 



Note that 5(0, z, x) = z. The function F q / g is independent of y and x, while F q j q has a 
simple, x independent, limit as y — ► 0. We have 

F, /9 (0,z,a;)=z 2 + (l- Z ) 2 , 

1 + z 2 ( 3 - 22 ) 

F, /9 (0,Z,X) = — - . 

Note that these limiting forms are the DGLAP kernels. In fact, one might say that the 
limit of the splitting functions as y — > at fixed z defines the DGLAP kernels for final 
state evolution. 

Using this result, performing the trivial ^-integration, and taking the y — > limit, we 
have 

(xJ a \Hl(t) - V(t)\{p,f,s',C,S,c} m ) ~ 

2pi ■ Q 



E 



2vr 



dz / dx 6 x 



Q 2 



({s'} m \{s}m) ({c'} m |{c} r 



1 



X 0(/, = g) T R [^ + (l-2)^] -5 X 



(3.23) 



+ = /a) C F 



1 + z 2 
1 - z 



1 



«5 x 



<5(x — x) 



3.5 The DGLAP equation 

We insert Eq. (3.23) into Eq. (3.2) and use the definition (2.2) of (x,f\p{t)) on the right 
hand side. We obtain 



d_ 
dt 



(x, U\p{t)) « g jT'dz |t r [z 2 + (1 - z) 2 ] 1 (x/z,g|p(t)) 



+ C F 



1 + z 2 
1 - z 



(x/z,/ a |p(i)) - (x,/ a |p(t)) 



(3.24) 



This is typically written 
d 



dt 



/a|p(*)) 



a s <iz 
2vr J £ z 

+ C F 



Tr [z + (1 — z) ] (x/z,g\p(t)) 

1 + z 2 ' 
1 - z 



(x/z,f a \p(t)) 



(3.25) 



The variable z is called x in Ref. [13], Eq. (2.21). However, in this paper, we use x for a different 
purpose. 
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This is the DGLAP equation for (x, / a where / a is a quark flavor. It has one term for 
g — > g + q splitting and one for g — > g + g splitting. The "+" prescription corresponds to 
the subtraction Eq. (3.24). It comes from the virtual splitting operator V(t) in Eq. (3.1). 



3.6 Transverse momentum ordered shower 

We have chosen virtuality as the shower ordering variable. One could use transverse mo- 
mentum instead, 

T l ({pJ} m+1 ) = log (-— %—„ ) . (3.26) 

V 2-2(1 - Z)pi ■Pm+lJ 

In this case, it is 

y T = z(l - z)y = -£f— e"* (3.27) 
2pi ■ Q 

that is fixed by the delta function 6(t — T{). The strong ordering limit is then y? — > at 
fixed z. In this case, one should include a factor 

9(z(l - z )> yr/ymax) (3.28) 

that limits the z integral in both the real and virtual splitting terms. However, this factor 
disappears in the limit yT — ► 0, leaving us with the same result (3.24). 



4. Soft gluon emission 



In this section, we examine the effects of soft, wide angle gluon emission. In Sec. 3, we 
argued that soft gluon emission could be ignored. Thus we dropped the interference term 
containing wik and, for the direct term wu for a q — > q + g splitting, we took the limit y — > 
with fixed z. The argument was that for a soft splitting, the observed quark feels no recoil, 
so that the virtual splitting term cancels the real splitting term. However, Ref. [4] argues 
that a dipole shower could fail to reproduce the DGLAP equation when the quark recoil 
is accounted for. Accordingly, we examine the contribution from soft, wide angle gluon 
emission when we account for the recoil using the momentum mapping from Ref. [1]. We 
write the q — > q + g contribution exactly, then examine the contribution in the soft emission 
limit. We find that this contribution vanishes in the limit of strongly ordered virtualities, 

y -o. 

We begin with Eq. (3.6) for the matrix element of TCi(t) — V(t). There is a term for 
a g — ► q + q splitting and a term for a q — > q + g splitting. We examine the q — > q + g 
term. We use the definitions for d( p and the evolution variable t, Eq. (3.15), as before. 
We introduce an integral of a delta function that defines the variable x as before. We also 
introduce an integral of a delta function that defines the angle i? between pi and p&. This 
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gives 



(x, h\Hl{t) - V(t)\{p, f, S ', C', S, C} m ) q ^ q+g 

yO(y < 



-1 

/ dx 8 1 x 



2m -q 



/0 

JlJm + 1 

1 A 1 X ( Q 1 i PrPkQ 2 

d cos v o cos w-H — — 

-i V PrQ Pk-Q 



o v Q 

X #(// = /a, /« = /a, /m+1 = g) 

X ({c'} m \t k (fk ~> A + / m +l) tj(/i ~* /i + /m+l)|{c} m ) 
24({p}m+l)%({/)?Wl) - Wil({/)P}m+l) 

o x 



(4.1) 



x — <5(x — x) 
PrQ ) 

4.1 The splitting functions 

For Wii, we use the exact expression given by Eqs. (3.19) and (3.20). For wik, we define 



2 A ik wik 



(4.2) 



Pi-Qy 

The (/ — > q + g interference function F[ n t includes the function Aik that describes how we 
divide the l-k interference term between a contribution that we associate with parton I and 
another that we associate with parton k. There is some freedom in choosing Aik- For our 
purposes here, we choose the version specified in Eqs. (7.2) and (7.12) of Ref. [14]. (Any 
of the other choices based on Eq. (7.2) from Ref. [14] would do as well.) The function i*j n t 
is simple when written in terms of the momenta pi, p m +i, pk and Q. When expressed as a 
function of x, i? and the splitting variables y, z, <j) according to the momentum mapping of 
Ref. [1] , F; nt is a little complicated. We find 

F mt (y,z,(j);x,tf) = 

2 2[gi(y, 1 - z,x,$) + gi(y,z,x,$)]g 3 (y,z) + 2yg 2 (y,x,$) 

gi(y,z,x,-d)g 3 (y,z) + yg2{y,x,$) - 2 tan(i?/2) v / xz(l - z)yg 3 (y,z) coscj) ( 4 - 3 ) 
2z 2y 
+ 1^2 ~ x{l- 5) 2 (1 + y) 2 ' 



where z was defined in Eq. (3.21), and where 
gi(y,z,x,tf) 



2^^y [1 - V ~ A] + \x tan 2 (tf/2) (1 - z) [1 - y + A] , 



g 2 (y,x,#)= 2 {1 + V)X , XX + \ tan a (tf/2)[2-(l + y)» + Ax] , (4-4) 
2(1 — x) 2 

g 3 {y, z) = Z - [1 + y + A] + [1 + y - A] . 

A simple calculation from Eqs. (4.3) and (4.4) shows that F[ nt is simple in the collinear 
limit, y — > with z fixed: 

F bA (0,z,<f>;x,&) = . (4.5) 
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This, of course, was obvious from the beginning. In a physical gauge like the Coulomb 
gauge used for our splitting functions, the interference graphs do not give a leading collinear 
singularity. Because of Eq. (4.5), we neglected the Wik term in the previous section. In this 
section, we do not simply take the limit y — > limit with z fixed. 



4.2 Kinematics 

The delta function 5(x — (pi • Q/pi • Q) x) in Eq. (4.1) relates x = 2pi-Q/Q 2 to x = 

%Pl ' Q/Q 2 an d the splitting variables y and z according to Eq. (3.16). We can solve 
this relation for x, yielding 



x = X(y,z,x) 



(4.6) 



The exact relation is rather complicated, 



X(y,z,x) 



2z(l-z)(l + yf 



(l + y)x-y(2z- if 



(2z - 1) + y)x - y] 2 - 4y 2 z(l - z) 1 . 



(4.7) 



As y — ► 0, X(y,z,x) approaches a very simple result X(0,z,x) = x/z. We rewrite the 
delta function as 



6(x - (pi ■ Q/pi -Q)x) = dX{y di Z,X) S(x - X(y, z, x)) 



(4.8) 



4.3 Evolution equation 

Using Eq. (4.8) and the substitutions for wik and wu, Eq. (4.1) becomes 

(x, h\Hi{t) - V(t) | {p, /, s>, c',s, c] m ) q ^ q+g = 

12 S J 0(yV < ymax ' ) A / dz ftn (^ s '^ m \^ m ) 
2pi-Q 



x J dx 



1 a i x ( q 1 , PrPkQ 2 
a cos it o cos v — 1 H — — 

i V PrQPk-Q 

t/ 



X e(fl = f a )(W}m\tk(fk -> /fc + g) *, T (/l /l +g)|{c}m> 

^int(y, 2, </>; sc> ^) - F q/q(V, z > x) 
dX(y,z,x 



dx 



5(x — X(y, z, x)) — 5(x — x) 



(4.9) 



We insert this into Eq. (3.2), giving the q — > q + g contribution to the evolution of the 
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energy distribution function. This contribution then takes the form 



d_ 



dz 



o 

x — A 
2vr 



2vr 



cZcos <d 6(y < y n 



Fmt{y, z i 4>; x, •&) - F q / q (y, z, x) 
9X % Z,£) (Xty,z,x),fl,#\p(t)) 



(4.10) 



where (x, f a ,mp(t)^ is a two parton correlation function, 
(x,f 3j ,'d\p(t)) = X^^i/ [ d {P>f' s '> c '> s > c }™] 



l k+l 



2pi-Q \ J , Q 2 Pi-Pk 

x — I o cos v — 1 + 



(4.11) 



Q 2 J V Q-pi Q-Pk 

X B(fl = U) ({s'}m\{s} m ) 

X ({<0m| W* - /* + g) 4 (/» /« + g)K C i™> 

x ({p,/,s',c',s,c} m |p(t)) . 

Thus we find an evolution equation similar to Eq. (3.24) except that the evolution of the 
one parton distribution is expressed in terms of a more complicated function, a two parton 
distribution. 

4.4 The soft emission limit 

If we were to follow the analysis of the Sec. 3, we would take the limit y — > at fixed 
z in the right hand side of Eq. (4.10). However, we should be a little careful about this. 
The kernel is singular not only for y — > at fixed z but also for y — > with (1 — z) = £y 
with £ fixed. This is the limit of soft, wide angle gluon emission when expressed using the 
variables y, z. Define the gluon emission angle i9 m +i by 



1 - cos •dm+x 



Then £ is related to $ m +i by 



1 - z 
y 



l + y 2 + (l + y)A 



2y 



Pm+l ■ Pi Q 2 

Pm+i ■ Q pr Q 



tan 2 (tf m+1 /2) + 



2y 



1 + cosi9 m+ i 



-i 



(4.12) 



(4.13) 



The range of £ is < £ < 1/y, with £ — > when i9 m +i — ► n and £ — ► 1/y when # m +i — ► 0. 
For y — ► at fixed i? m +i, 



1 



(4.14) 



tan 2 (i9 m+1 /2) 

Our purpose in this section is to investigate whether the soft emission integration region, 
y — > at finite £, contributes to the evolution at small y. 
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The integration over the soft region can be expressed by using £ as the integration 
variable, with 

dz = yd£ . (4.15) 
Taking y — > at finite £, the factor y from the integration measure times the kernel becomes 



F int (y, z, 4>; x, i?) - F q/q (y, z, x) 



2xtan 2 (i9/2) 



(4.16) 



2 + (l + x£)tan 2 (^/2) - 2 tan(??/2)^cos 



Thus we would have a finite contribution to the integration from the soft emission region 
in the limit y — * if we were to multiply the splitting function by anything that is finite 
in the soft emission limit. 

However, for our inclusive observable, we need to multiply by the function 



dX(y,z,x) 



dx 



(X(y,z,x),L$\p(t)) - (xj a ,#\p(t)) 



Using the expansion 



(4.17) 



X{y,z,x) ~ x + y(l-x + tx) + 0(y 2 ) , 



(4.18) 



we see that the factor (4.17) becomes in the soft limit 



d_ 
dx 



(l-x + £r)(x,/a,4o(t))l +0(y 2 ) 



(4.19) 



Thus the integrand has a net factor y in the soft limit, so that the soft region does not 
contribute to the leading shower evolution of (x, / a (We have investigated the q — > 

q + g contribution here, but the g — > q + q contribution vanishes in the soft limit.) 

We note that the fact that the soft region does not contribute follows from the inclusive 
nature of our observable. With other observables, there would be a leading soft gluon 
contribution. 

We note also that there is an angular ordering cancellation contained the soft limit 
splitting kernel, Eq. (4.16). If we did not have the factor (4.17), the important integration 
region for £ when tan 2 (i?/2) <C 1 would be £tan 2 ($/2) > 1- That is, tan 2 (tf m+ i/2) < 
tan 2 ($/2). However, this effective angular ordering limit on the £ integration does not play 
a role in the present analysis. 

This analysis has covered a virtuality ordered shower. For a transverse momentum 
ordered shower, the variable yx = z(l — z)y that was introduced in Sec. 3.6 is fixed by the 
evolution time. The wide-angle soft limit is yx — > at fixed £' where (1 — z) = yJWyr- 
Then the analysis just given for virtuality ordering applies also to transverse momentum 
ordering with ■Jyr playing the role of y. 
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5. Catani-Seymour shower scheme 



We now examine what happens with a Catani-Seymour dipole shower as defined in Ref. [9] . 
A parton with momentum pi splits to daughter partons with momenta pi and p m +l- There 
is a helper parton with momentum p^ that absorbs some momentum from the splitting, so 
that after the splitting its momentum is p/.. The splitting is described by splitting variables 
y, z, and <p. In contrast to the momentum mapping of Ref. [1], in the Catani-Seymour 
scheme, the momenta of all other partons are unchanged by the splitting, pi = pi for 
i $l {I, k}. The momenta pi, p m +i, and pk are given in terms of pi, and pk and the splitting 
variables by 

Pl = zpi + y(l - z)p k + [2yz{\ - z)pi -p fc ] 1/2 KLL , 
p m+1 = (1 - z)pi + yzp k -[2yz(l - z)pi-p k ] 1/2 K± , (5.1) 

Pk = (i - y)Pk ■ 

Here k± is a spacelike unit vector orthogonal to pi and p k that has azimuthal angle <j). Note 
that y is proportional to the virtuality of the splitting 

Pi' Pm+l /- s 

y = . (5.2 

Pi - Pk 

For the observables that we need to analyze, we need some additional variables involv- 
ing Q, 

a _ 2prQ 

X ~ Q2 ' 

_ ggrQ 
X ~ O 2 ' 

n2 ( 5 - 3 ) 
pi • Pk Q 

PrQ Pk-Q 
Pk ■ Q 

r = . 

Pi ■ Q 

Using the first equation in Eq. (5.1), we find how x is determined by x, r, and the 
splitting variables, 

x = Z(z,cp,y;r,i}) x , (5-4) 

where 

Z(z, 4>, y,r,$) = z + y(l — z)r — \J yz(l — z)r sin($/2) cos <fi . (5-5) 

Note that Z(z, <fi, y, r, ■&) — > z for y — > 0. 

In defining the shower evolution, we take a spin average everywhere, so that spin 
indices disappear. In addition, we use the leading color approximation. Then we have just 
one color label for the m parton state, {c} m . The evolution equation is 

d 



1 — COS "& 



— [X 



fMt))=/Z^ [ MP, Mm] 

m 

x(x,f a \H? s (t)-V cs (t)\{p,f,c} m ) (5 - 6) 

X ({Pj,c}m\p(t)) , 
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where we integrate over a complete set of parton states \{p, f, c} m ) just before the splitting 
represented by the real and virtual splitting operators Tif (t) and V cs (i). 

The matrix element of the splitting operator for a Catani-Seymour shower has the 
form 



{xJ a \H? s (t) - V cs (t)\{pJ,c} m ) 



2pi ■ Pk 

167T 2 



EE E 

■m+l 



dy (i - y) / dz 



2- 



'/.,/< " ^ Q2 

L J_ 

x Xik{{c} m ) wg s ({f,p} m+1 ) . 



2tt 



2pi-Q 



(5.7) 



For each splitting of a parton I we choose a partner parton k. The momentum mapping 
depends on both I and k, so the sum over k comes to the outside. There is a sum over flavor 
choices Q for the splitting. If fi is a quark or antiquark flavor, there is only one choice, 
fl —* fl + g- When fi=g then the gluon can split to two gluons or to any quark-antiquark 
pair. There is an integration over the splitting variables y, z, 4> with the jacobian that 
follows from the Catani-Seymour kinematics. There is a delta function that defines the 
shower time t. We will make the simplest choice of a virtuality ordered shower, 



mpj} 



m+l) 



log 



Q 2 



2pl ■ Pm+l 



log 



Q 2 



y %pi ■ Pk 



(5.8) 



Thus the y-integration can be eliminated, with 



Ql 



%Pl ■ Pk 



(5.9) 



In the next factor in Eq. (5.7), the first term is from the splitting operator, 7i^ s (t), while 
the second term is from the virtual splitting operator, V cs (t). Both operators contain the 
same splitting function, represented in the last line. 

In all cases except a g — > q + q splitting, the splitting function includes the possibility 
of emitting a gluon with label m + l from a parton with label % and emitting the gluon from 
a parton with label j, treated coherently. The total emission probability is divided into a 
term associated with parton i, for which we use / = i and k = j and another contribution 
associated with parton j for which we use I = j and k = i. The splitting function is w^jp. 
For q —¥ q + g, 



™?k S 



Aixa s 

ypi - Pk 



C F 



1 — z + zy 



H + z) 



(5.10) 



Finally, Xik({c} 

Xlk({c}m) = 0. 



1 if partons / and k are color connected in color state {c} m . Otherwise 
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We now specialize to the case that / a is a quark flavor. The case of a g — ► q+q splitting 
is not particularly instructive, so we consider only ag->g + g splitting. Then 



(*,/a|p(*)) 



Jo 2vr 2vr v y; [i- z + zy 

dcosfi / dr 

(x/Z(z, 4>, y; r, / a , r, 
Z(z,4>,y;r, •&) 



- (xj a ,r,#\p(t)) 



(5.11) 



where 



x (5 x - 



<5 2 

x ({p,/,c} m |p(*)) . 



5 COS 1? - 1 + 



Q 2 gj • Pfc 
<5 • pi Q-Pk 



5 r- 



Pk-Q \ (5-12) 

pi ■ Q 



Here (a;, / a , r, $|p(t)) is a two parton correlation function. One parton, with flavor / a , 
carries energy fraction x. Another parton, color connected to the first, is separated from 
it by an angle i? and carries energy fraction rx. 

We can now make the approximation of strongly ordered virtualities in the shower. 
For this, we simply take the y — > limit. In this limit, the splitting function simplifies and 
Z(z, 4>, y; r, #) ~ z. Then 



d_ 
dt 



(x,h\p(t)) 



2tt 



d(f> a s 



dz — C F 

/o 2vr 2vr 



1 - z 



J dcosd J 



dr 



l/zj & ,r,0\p(t)) - (x,h,r,#\p(t)) 



(5.13) 



Now nothing depends on <fi, so we can trivially perform the (^-integration. Furthermore, 
the integral over cosi? and r of the two parton correlation function is the one parton 
distribution function, 

/l poo 
dcostf J dr {x,U,r,0\p(t)) = (x,f a \p(t)) . (5.14) 

Thus we obtain, within the strongly ordered virtualities approximation, the q — > q + g 
contribution to the DGLAP equation for the energy distribution function of a quark, 



(a,A|p(t)) 



Qs 

2tt 



1 1 + z 2 

dz Cf 



I q^q+g " n J 

This is often written using the notation 



l-z 



1 



(x/zj a \p(t))-(x,f a \p(t)) 



(5.15) 



d 
dt 



%U\p{t)) 



dz 



2tt Jo z 



C f 



1 + z 2 
1 - z 



(x/zj a \p(t)) . (5.16) 
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We have presented the derivation for a virtuality ordered shower. The Catani-Seymour 
dipole showers of Refs. [6, 7, 10] use transverse momentum as the ordering variable, 



Ti({p,f} 



m+l) 



log 



Q 2 



Then the variable 



y T = z(l - z)y 



2z{l - z)pi -pm+1 

Q 2 



(5.17) 



(5.18) 



2pi • Pk 

is fixed by the delta function <5(i — TJ). The strong ordering limit is then yx — > at fixed 
z. In this case, since the definition of y entails that y < 1, one should include a factor 



>(1 - z) > y T ) 



(5.19) 



in the integrand. However, this factor disappears in the limit yx — > 0, leaving us with the 
same result (5.15). 

In our derivation, we have simply taken the limit y — > (or yx — ► 0) inside of the 
integration over z. One could explicitly examine the region in the z-integration that cor- 
responds to soft gluon emission, y — > with (1 — z)/y = £ finite. To do that, one changes 
variables from z to £ and takes the y — ► limit inside the integration over £. This leads to 
what would be a finite contribution to the integral except for the factor 



(x/Z(z, 0, y; r, / a , r, 
Z(z,(/),y;r, ti) 



{xj a ,r,0\p(t)) 



(5.20) 



For a virtuality ordered shower, this factor provides a factor y that removes the soft emission 
contribution in the y — > limit at fixed £ = (1 — z)/y. For a transverse momentum ordered 
shower, this factor provides a factor ^fyr that removes the soft emission contribution in 
the j/t — ► limit at fixed £' = £/z = (1 — z) 2 /yT- This justifies taking the limit y — > (or 
— * 0) inside of the original integration over z. 



6. Conclusions 



In a final state parton shower, the distribution of the energy carried by a single parton 
evolves as the shower evolves and the energy carried by each mother parton is distributed 
among the daughters. We have examined this evolution for a final state parton shower 
in electron-positron annihilation. Our principle example is based on the shower evolution 
equations of Ref. [1] , which include complete color and spin information. These evolution 
equations are based on a dipole picture of parton splitting, in which gluon emission from a 
parton I can interfere with gluon emission from another parton k. We have also analyzed 
evolution in (spin averaged, leading color) shower evolution based on the Catani-Seymour 
splitting functions and momentum mapping, as used in the practical shower programs of 
Refs. [6, 7] and in the final state showers of Pythia [10, 11]. 

The use of a parton shower to model the evolution of the final state created from a 
short distance scattering process is based on the idea that successive steps in the evolution 
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are separated by happening on very different time scales. That is, each successive step 
involves much smaller virtuality (or transverse momentum) than the previous step. The 
shower evolution is based on factoring the amplitude for a relatively soft splitting from the 
previous harder splittings. The actual evolution defined by the evolution equation includes 
integrations over the virtuality of each splitting, including an integration region in which 
the virtuality is not much smaller than the previous virtuality. 6 Whatever is left over 
coming from the region in which the virtualities of successive splittings of the same parton 
are similar rather than strongly ordered may be considered to contribute to the second 
order DGLAP kernel. Thus this contribution is suppressed by a factor of a s . 

We have examined the evolution of the distribution of the energy carried by a parton, 
[x, / a \p{fj) , as defined in Eq. (2.2). The analysis presented in Sec. 3 for the shower evolution 
equations of Ref. [1] covered the evolution of the energy distribution for quarks. For the 
sake of completeness, we include the evolution of the energy distribution for gluons in 
Appendix A. For evolution using the Catani-Seymour splitting functions and momentum 
mapping, we have limited the analysis to the evolution of the energy distribution for quarks 
from q — > q + g splitting. In each case, we emphasized a virtuality ordered shower, but 
pointed out that a transverse momentum ordered shower gives the same result for the 
evolution of the parton energy distribution. In each case examined, we find that the energy 
distribution obeys the standard lowest order DGLAP equation when we approximate the 
evolution using the strong ordering limit. This result addresses the concern expressed by 
Dokshitzer and Marchesini [4] that evolution in dipole based showers might fail to follow 
the DGLAP equation. 

The shower evolution time t represents the scale fi 2 = 2pi • p m +i at which one parton 
is resolved into two daughter partons. When /j, 2 is much smaller than the scale Q 2 of the 
hard process, the logarithm t = log(Q 2 /// 2 ) is large. When (x,/ a |/9(t)) is expanded in 
powers of a s , the perturbative coefficients contain powers of t. The evolution equation 
enables us to sum the large logarithms of Q 2 / /U 2 . There are a number of other cases in 
which a physical observable presents us with a large ratio of physical scales and in which 
it is known how to sum the large logarithms of this ratio. In some of these cases, parton 
shower evolution will generate the correct summation. This is what we found here for the 
parton energy distribution, which, although not directly observable, maps to the observable 
hadron energy distribution. In some other cases, it is pretty straightforward to see that a 
parton shower generates the correct summation, as for the dependence of jet cross sections 
on a jet resolution parameter based on kx when the jets are generated with a %• ordered 
parton shower. In other cases, parton shower evolution will generate a summation of large 
logarithms, but it may be the wrong summation. For instance, we expect that standard 
parton showers do not generate "BFKL" logarithms that describe the physics of large 
rapidity differences, as in [15]. 

For a given choice of shower evolution it is of significant interest to know what sum- 
mations of large logarithms the shower evolution equation correctly produces. That is, 

For this reason, one attempts to define the splitting functions and momentum mapping in such a way 
that the parton shower model does something sensible when the ratio of successive splitting virtualities is 
not small. 
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one would like to validate particular shower evolution schemes against known results for 
summations of large logarithms. Consider, for example the energy-energy correlation func- 
tion in electron-positron annihilation at nearly back-to-back angles 6. One knows the 
summation of large logarithms log(l/^ 2 ) in full QCD. Then one needs to derive the cor- 
responding summation in the shower model, deriving the appropriate evolution equation 
from the shower evolution equation. A second example in electron-positron annihilation is 
the summation of logarithms of 1 — x in the energy distribution of partons or hadrons. For 
hadron-hadron collisions, one should investigate, for example, the transverse momentum 
distribution of Drell-Yan lepton pairs. 

We note that the existing literature goes a long way toward validating parton shower 
evolution against known large logarithm summations in the case of parton showers based 
on using the angle between daughter partons in a splitting as the shower evolution variable. 
Indeed, some of the derivations of large logarithm summations in QCD make use of on an 
angle-ordered parton cascade picture [16, 17, 18, 19, 20]. Correspondingly, the paper [21] 
argues that a parton shower ordered in splitting angles gets the leading terms in certain 
large logarithm summations correctly. 

To our knowledge, less is known analytically about how well virtuality or kx ordered 
showers reproduce known large logarithm summations. We do note that Ref. [22] looks at 
large logarithm summations for both angle ordered and for virtuality ordered or kx ordered 
showers with respect to "non-global" observables. 

To build from the existing results, one would need to check carefully, for a given shower 
scheme, the accuracy (in powers of a s and the logarithms) with which the full QCD result 
for a given kind of summation is reproduced by the shower. As emphasized by Dokshitzer 
and Marchesini [4] , one would need to pay close attention to the treatment of momentum 
conservation in the shower. 

We do not attempt a general analysis here. What we offer in this paper applies only 
to the circumscribed problem suggested in Ref. [4]. For further investigations along these 
lines, we believe that it is useful to have shower evolution expressed as a definite integral 
equation, as in Ref. [1]. We hope that the present paper provides some hints about how a 
program of validating parton shower evolution against known large logarithm summations 
could be pursued. 
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A. Energy distribution of gluons 

We have analyzed the evolution of the energy distribution (x, / a Lo(i)) in the case that / a is 
a quark flavor. The case of an antiquark flavor is so similar that it hardly needs a separate 
treatment. However, there are some differences in the case that / a = g. We provide a brief 
analysis of this case in this appendix. 

There are three sorts of contributions. The first is a q — ► q+g or q — ► q+g splitting, as 
depicted in Fig. 6, in which it is the energy of the daughter gluon, instead of the daughter 
quark or antiquark, that is measured. The second is a g — > g + g splitting, as in Fig. 4 
but with the q + q replaced by g + g. Either of the two daughter gluons can be the one 
whose energy is measured. There is a virtual splitting counter term for this contribution. 
There are g — > g + g interference graphs, but these are singular only in the soft limit and 
do not give leading contributions for the evolution of the parton energy distributions, as 
we have seen in some detail for quark splitting. (Evidently, if it is the soft gluon itself that 
is measured then there is no contribution to the evolution of (x, f & \p{t)) at any finite x 
in the limit in which the soft gluon energy vanishes.) The final contribution is the virtual 
splitting counter term for a g — > q + q splitting, as depicted in Fig. 4. 

For the gluon case, Eq. (3.18) becomes 

(x,g\Hi(t)-V(t)\{ P ,f,s',c',s,c} m ) « 

if 8^ H i 2«L \ ) 

fljm+l 

X ({s'}m|Mm} ({ c '}m\{c}m) 

x M/i ¥=s,fi = fiJm+i =g) c F wu({f,p} m+1 ) y~, s ( x ~ rzz) ( A - 1 ) 

+ Q(fl =gjl = gjm+l = g) C A W U ({f,p} m+1 ) 

X — - — S\X — ) H S\ X — — ) — S(x — x) 

1 — z \ 1 — z ) z \ z ) 

- 0(fl = g,fl € Q, fm+1 = ~fl) T R Wu({f,p} m +l) S(x - X) 

There are five terms in three groups. In each case, we have adopted the notation that 
parton m + 1 carries momentum fraction 1 — z. In the first term, the mother parton with 
flavor fi is a quark or antiquark that splits into a parton of the same flavor with index I and 
a gluon with index m + 1. We measure the energy of the daughter gluon, so x = (1 — z)x. 
In the next three terms, the mother parton I is a gluon and daughter partons I and m + 1 
are both gluons. The definitions of Ref. [1] explicitly break the symmetry between gluons 
I and m + 1 in such a way that there is a singularity in wu corresponding to gluon m + 1 
being soft, but no singularity corresponding to gluon I being soft. In the second term, we 
measure the energy of the daughter with label m+ 1, so x = (1 — z)x. In the third term, we 
measure the energy of the daughter with label /, so x = zx. The fourth term is the virtual 
splitting counter term, with x = x. The fifth term is the virtual splitting counter term for 
a g — ► Q + Q splitting. Here fi £ Q means that fi is a quark flavor, Q = {u, d, s, . . . }. 
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All that we need is the splitting functions wu in the y —* limit. For fi^gwe have 

_ Aira s 1 + z 2 



Pi-Q y 1 - z 



(A.2) 



as in Eqs. (3.19) and (3.22). For ft = g, fi = g, we have, from Eqs. (2.46) and (2.52) of 
Ref. [13] 7 , 



4-7ra s 

prQ y 



2z 
1 - z 



+ z(l-z) 



For fi = g, fi e Q, we have 



wu 



47TQ S 

vi- Qv 



[z 2 + {l-zf] , 



(A.3) 



(A.4) 



as in Eqs. (3.19) and (3.22). We make these substitutions and perform the trivial <j) integra- 
tion. In the first two of the five terms, we use the freedom to change integration variables 
to interchange the roles of z and 1 — 2. This gives 



(x,g\Hi{t) - V{t)\{p, /, s 1 , c, s, c} m ) 

2prQ 



y / dz dx 5 
z Jo Jo 
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(A.5) 



5(x — x) 



In the last line, rif is the number of quark flavors, that is, the number of elements of Q. 

We can now substitute this result into Eq. (3.2) for the evolution of (£, g|/o(£)). On 
the right hand side of the evolution equation, we use the definition (2.2) of (x, f a \p(t)). 
This gives 



d_ 



l + {l-z) 2 1 



[x/z,f'\p{t)) 



+ C A 
+ C A 



2(1-2) 



+ 2(1-2) 



22 
1 - 2 



+ 2(1-2) 



'- (x/z,g\p(t)) 
x/z,g\p(t)) - (x,g\p(t)) 



(A.6) 



T R n f [2 2 + (l-2) 2 ] (x,g|p(t))| 



7 We interchange z and (1 — z) to match the conventions of this paper. 
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This is the DGLAP equation. It is usually written as 



d_ 
dt 




l + {l-z) 2 



^ + <i-.>) 



z 



(x/z,f'\p(t)) 



(A.7) 



llC A -4r R n f 
6 



6(1 -z) (x/z,g\p(t)) 
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